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1. Introduction 

Theoretical arguments and experimental observations indicate that new particles or inter- 
actions play an important role at the TeV scale, which will become directly accessible at the 
Large Hadron Collider (LHC) - scheduled to start in 2008 - and its planned complement, 
the International Linear Collider. In the near future we can therefore anticipate ground- 
breaking discoveries that reveal physics beyond the Standard Model (BSM) and allow to 
gain insight into the structure of the fundamental theory. Theoretically appealing exten- 
sions of the Standard Model (SM) often feature numerous additional interacting heavy 
particles. Supersymmetric theories [1], for example, are attractive, because they solve the 
hierarchy problem and allow for the unification of electroweak and strong interactions. The 
Minimal Supersymmetric Standard Model (MSSM) is one of the best studied candidates 
for BSM physics. Its phenomenology is characterized by sparticle production and cascade 
decays, which lead to many-particle final states and scattering amplitudes with complex 
resonance structure. Cascade decays also occur in other extensions, e.g. in universal extra 
dimensions models [2]. 

In order to extract the additional Lagrangian parameters of an extended theory from 
collider data, theoretical predictions are required that match the experimental accuracies. 
This can usually only be achieved by taking into account higher order corrections in per- 
turbative calculations. Next-to-leading order calculations for phenomenologically relevant 
2 — ► n processes with n > 4 are technically very challenging or not yet feasible [3, Sec. 30]. 
Consequently, production and decay stages are regularly factorized by means of the narrow- 
width approximation (NWA), which effectively results in on-shell intermediate states. 1 Its 
main advantage is that sub- and nonresonant as well as nonfactorizable amplitude contribu- 
tions can be neglected in a theoretically consistent way. Huge calculational simplifications 
occur already at tree level. For these reasons, the NWA is employed in nearly all studies 
of BSM physics. Note that it is implicitly applied whenever branching ratios are extracted 

lr rhe NWA can thus not be applied if on-shell states are kinematically forbidden. 
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from scattering cross sections. A reliable NWA uncertainty determination is therefore cru- 
cial. Given the width T and mass M of an unstable particle, the uncertainty of the NWA is 
commonly estimated as 0(T/M) for each Breit-Wigner propagator that is integrated out, 
with r/M typically < 5%. For larger widths nonresonant contributions can no longer be 
neglected. 

Recently, two circumstances have been observed in which the standard NWA is not 
reliable [4,5]: the first involves decays where a daughter mass m approaches the parent mass 
M; the second involves the convolution of parton distribution functions with a resonant 
hard scattering process at center-of-mass energy y/§. In this article we elucidate that both 
effects arise due to a significant deformation of the Breit-Wigner shape that is caused by 
threshold factors, and is not restricted to the region where the Breit-Wigner is cut off, 
i.e. where M — m or \/l — M is approximately T. An essential factor is that the amplitude 
can contribute additional powers of the threshold factors, which strongly amplifies the 
effects. For sample applications we then demonstrate that 0(r/M)-accurate predictions 
can nevertheless be obtained by integrating out the Breit-Wigner in combination with the 
relevant threshold factors. 



2. NWA modifications 

To illustrate why the NWA error becomes unexpectedly large for mass configurations in 
an extended vicinity of kinematical bounds and how it can be modified in such cases, we 
consider the partial decay rate of a heavy particle A that predominantly decays in two 

1 2 

stages via an intermediate resonance C, i.e. A — ► B,C and C — > D,E. In terms of the 
n-body phase space element 

n n ,3 

#(P; P1 , ... tPn ) = (2^)V 4 )(P - 5» ft TT^TW ^ 

i=i i= i [A7T > 

and the matrix element M, the off-shell decay rate is given by 

r o fr-shcii = / #l-^! 2 (2-2) 

= 2M~ A J § D{p2c) J dMp2c) J dMp2c) l-^^l 2 • ( 2 - 3 ) 

In Eq. (2.3), the phase space factorization 

d(f) = dfa^dfo (2.4) 

has been applied, where d<pi{d(j)2) is the 2-body phase space element of the first (second) 
decay stage. In the rest frame of A, 

#1 = lk* '^t ^ ^ \Pc\ = ^Y^ M B + \/pi,M A )P(M B -^p%,M A ) , (2.5) 
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where (3(m, M) = \J\ — m 2 /M 2 . For d(f>2 one finds an analogous expression. In addition 
to d(f>, also \M.\ 2 has been factorized into the squared propagator denominator 

s (rt^wrwi (2 ' 6) 

with 4-momentum pc and |A4 r | 2 , the residual squared amplitude for the A —> B,D,E 
decay. Starting in Eq. (2.3), we have indicated the particularly relevant ^-dependence of 
quantities explicitly, but suppressed the dependence on other kinematical variables. In the 
limit — > 0, D{pq) is asymptotically equal to 2irK5(p 2 c — M c ) with 

K = = I" ^D(q 2 ) . (2.7) 

2M C T C J-oo 2vr W ' K ' 

This replacement constitutes the standard NWA. Employing it, one obtains 

Tnwa = ^-K J #i(M 2 ) J d(j> 2 (M c ) \M r {M 2 )\ 2 . (2.8) 

It is suggestive to mitigate threshold-induced deviations by absorbing the amplifying factors 
of /3-form that occur in d<j>i (see Eq. (2.5)) and d(f>2 into K. The absorbed factors have to 
be normalized to p 2 -, = Mq, so that they are not taken into account more than once. K is 
thus replaced by 



»(M A -M B ) 2 dp 2 
H Me? 271 



J(M D - 



P(M B + ^p 2 c ,M A )(3(M B - ^p 2 c ,M A ) 

(3(M B + M C , M A )(3{M B - M C , M A ) 
P(M D + M E , yfpl)P(M D - M E , yjtf,) 

(5{M D + M E , M C )P(M D - M E , M c ) 
f\Mr\<\Jph> M A,M B ,M D ,M E ) 
f lMrl 2(M c ,M A ,M B ,M D ,M E ) ' 



(2.9) 



Below we find that additional amplifying factors like M 2 — m 2 = (3 2 (m, M)M 2 can arise 
due to momentum-dependent residual matrix elements. Such factors are included gener- 
ically in Eq. (2.9) as f\M r \ 2 - They are process specific and their powers depend on the 
spin/polarization of external states. 

Note that replacing K with K does not affect the intrinsic properties of the NWA 
and that it generalizes to multi-body decays. A closed-form result for R = K/K can only 
be given for special cases. In Fig. 1 we show for the ratio R2 that is obtained by only 
taking into account the /3-factors that arise from dfo and setting the upper integration 
boundary to infinity the deviation from 1 normalized to Tc/Mq- One can see that the 
largest deviation occurs for + M B — ► Mq and M B (or Md) — > 0. The sizable 
effect for m = Md ~ M = Mc and a small mass M B , which we set to zero to obtain 
analytical results, is further amplified if the matrix element of the second decay contributes 
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Figure 1: The relative deviation of the modified from the standard NWA factor, normalized to 
Tc/Mc — 1%, i.e. (i?2 — l)/{Fc/Mc), is shown as function of the masses of the second decay 
C — > D,E (see main text for details). Contour lines are shown for the values (solid), —1 (dot- 
dashed) and 1, 3, 10, 50, 100 (dashed). The dash length decreases with increasing magnitude. 



an additional factor M 2 — m? . This is for example the case if C and E are fermions and D 
is a scalar (and spin correlations between decay 1 and 2 are neglected) or if C is a scalar 
and D and E are fermions. For these decay types, strong effects have been observed in 
Ref. [4] and Ref. [5], respectively. The corresponding R' 2 = K' 2 /K is given by 




1 



(q 2 — m ) /q 



2vr {q 2 - M 2 ) 2 + (MT) 2 (M 2 - m 2 ) 2 /M 2 

dq 2 1 

(q 2 - M 2 ) 2 + (MT) 2 

with T = Tc- In Fig. 2 we show the deformation of the Breit-Wigner shape due 
additional threshold factors when m approaches M. After integration, we obtain 



(2.10) 
to the 
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(2.11) 



with 7 = r/M, (3 = /3(m, M) and A = g^ ax /M 2 — 1 when factors of 1 +7 2 are approximated 
by 1. (This approximation does not produce a visible difference in Fig. 3.) The result 
confirms that away from threshold, where j3 ~ 1, one obtains R' 2 ~ 1 with 7 <C /? 2 ,A. 
When approaching the threshold, i.e. (3 — > 0, the divergence of the second, formally 7- 
suppressed term overcompensates the decrease of the first term. In Fig. 3 we show the (3 
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Figure 2: The Breit-Wigner shape deformation is displayed that is caused by threshold factors 
when a decay daughter mass m approaches the parent mass M. More specifically, the integrand of 
the numerator (solid) and denominator (dashed) of Eq. (2.10) are shown in unspecified normaliza- 
tion as functions of the invariant mass \fq*. T/M — 1% and m = M — 2T. 



dependence of R' 2 for typical values of \J q^ &x /M. The deviation from the standard NWA 
clearly exceeds 0(T/M) already for threshold masses m that are still significantly below 
the resonant region roughly bounded by M ± T. We note that if \J q^^/M — 1 < 7 
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Figure 3: The decay mass dependence of the relative deviation of the modified from the standard 
NWA factor is shown in units of the conventionally expected uncertainty T/M = 1% for yj q^ ax /M £ 
{1.05, 1.1, 2, 10}. The dash length decreases with increasing y^max- 
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the threshold amplification is confined to the resonant region. However, in this case the 
arctan terms are no longer approximately 7r/2, which results in a much larger than expected 
uncertainty of the standard NWA for arbitrary values of m. If g^ax ^ M 2 the contribution 
from the region q 2 q^&x t° R'2 ls enhanced by the factor (q 2 — m 2 ) 2 /(q 2 m 2 ) « q 2 /m 2 and 
the production cross section's suppression close to threshold can become important. 



3. Applications 

In this section we demonstrate that the NWA modification proposed in Sec. 2 allows to re- 
duce the uncertainty to the conventional expectation for mass configurations in the vicinity 
of kinematical bounds. 

As a first application we study the threshold-improved approximation for scalar scat- 
tering and decay. More specifically, we study the processes displayed in Fig. 4. For this 
type of process large standard NWA deviations have been observed in Refs. [4, 5] when 
m d ~ M or the center of mass energy y/s ~ M. We assess the quality of the modified 




Figure 4: Process 1 (left) with scalar scattering and decay and process 2 (right) with decay into 
non-scalar particles (fermions). Lines without labels correspond to massless particles. 

NWA by comparing the off-shell cross section to the cross section cttnwa calculated in NWA 
with K of Eq. (2.9), where M A = *fs,M c = M,T C = T,M D = m d and M B = M E = 0. 
The deviation is measured in units of T/M, which we set to 0.01, using 

^INWA - ( a ~^ ~ l) ■ (3.1) 

We start by neglecting matrix element effects and thus set f\M r \ 2 = 1- F° r process 1 with 
m p ~ M we find satisfactory NWA uncertainty reduction. For instance for m p = 1.1M, 
-Rinwa ^ 3 as long as (y^ — M)/M > 10~ 5 . For m p ^ M, however, large deviations 
occur in a significant parameter space region, in particular for ~ M. For process 2 
large deviations remain, independent of the value of m p . The Breit-Wigner deformation 
arises apparently not just from threshold-type phase space element factors. To achieve a 
satisfactory NWA uncertainty reduction it is in general essential to take into account factors 
originating from the matrix element that distort the Breit-Wigner shape. We separate them 
into production and decay-related factors: 

f\ Mr \2(M,^/s,0,m d ,0,m p ) = f p (M,y/s,m p ) f d (M,m d ) . (3.2) 
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For process 1, the decay matrix element is a coupling constant and we thus have f d = 1. 
For process 2, however, we have 

f d (M,m d ) = l -^f = ^-^=p\m d ,M)^ , (3.3) 
m d m d m d 

where we have divided by m 2 d to obtain a dimensionless quantity and expressed the result 
in terms of threshold /3-factors that also appear in the decay phase space element (see 
Eq. (2.5)) in order to show the deviation-amplifying effect of the decay matrix element by 
contributing an additional power to the Breit-Wigner deforming factor in d<p d . When taking 
decay matrix element effects into account, i.e. employing K with f\M r \ 2 = fd(M,m d ), 
NWA deviations are mitigated to 0(Y/M), except for the region yfs < 1.5M, where the 
^-channel production matrix element causes significant Breit-Wigner deformations. These 
production effects can be remedied with 

f p (M, y^, m p ) = (3- 2 y\M*-m% ■ (3.4) 

We now extend our analysis to more complex processes and study the NWA uncertainty 
reduction at hadron colliders for sparticle production and decay in the MSSM. In Ref. [4] 
the standard NWA accuracy was studied for the process ud —> (g —> !sl,rs)Xi ■ Here, 
the resonant particle, i.e. the gluino, is produced in a t-channel process with either d^ or 
ul exchange. For this process a variation of the sl mass between and the gluino mass 
revealed unexpectedly large NWA deviations for squark masses that are larger than 0.8M g . 
The slope of the increasing deviation when the squark mass approaches the gluino mass is 
qualitatively very similar to the slope displayed for R' 2 in Fig. 3. Since R' 2 does not take 
into account the t-channel production effects, we conclude that they do not significantly 
alter the dominant decay effects. We have confirmed that in this region the uncertainty 
of the NWA is reduced to 0(Y/M) if K' 2 of Eq. (2.10) is used with M = M g , Y = Y g , 
m = M-g and y^max = l-4Mg (matched at the squark mass value where 0" o fr_ s heii = <tnwa)- 
We note that for sl masses below 0.8M g the NWA overestimates the off-shell cross section 
by up to about 20%. This deviation is, however, consistent with an expected uncertainty of 
0(Y/M), since in this region the gluino width increases to about 10% of its mass. Ref. [4] 
also illustrates deviations for the sl~sr decay asymmetry, which are consistent with NWA 
corrections of 0(Y/M). 

As a last application we consider cascade decays, which are the natural testing ground 
for Eq. (2.9). More specifically, we study gv,L production at the LHC, i.e. in proton- 
proton collisions at 14 TeV, with the subsequent cascade decay g — > sls and sl — ► Xi c 
at the SPSla' benchmark point [6] in the MSSM parameter space. Phenomenologically, to 
consider a squark decay into the LSP candidate X\ would be more natural, but the resulting 
complete Feynman amplitude features a complicated resonance structure whose study we 
leave to future work. Even for the gluino decay chain considered here, interference arises 
from g — > (c* L — > X\^) c - Its effect is, however, small. We confirmed that omitting it does 
not affect our 0(Y/M) accuracy goal. In this article we focus on the resonant sl state 
(with M = 570 GeV and Y = 5.4 GeV at SPSla') and the NWA accuracy relative to Y/M 
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that is obtained using K with 

M? - Mf M| - M| „ M~ 

/ w (M ? ,M 5 ,0,M^,0) = g 5 2 X = P 2 (Ms,M~ g ) (3 2 {M^Ms) ^ (3.5) 

9 x x 

versus K when the strange squark mass approaches either the gluino or chargino mass of 
607 and 184 GeV, respectively. 2 Results calculated with MadEvent [7] and Sherpa [8] using 
CTEQ6L1 parton distribution functions [9] and spectra and decay widths obtained with 
SPheno [10] and SDECAY [11] are displayed in Fig. 5. The Monte Carlo integration error 
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Figure 5: The accuracy of the NWA cross section normalized to the conventionally expected 
uncertainty is shown for gv,L production at the LHC followed by the cascade decay g — > sls and 
sl — * Xi c m the MSSM at SPSla' for a variable strange squark mass that approaches the chargino 
mass (left) and the gluino mass (right). Results are displayed for the standard NWA (diamonds) 
and the improved NWA (INWA) of Eq. (2.9) (boxes). T(s L )/M(s L ) ranges from 0.03% to 0.16% 
(left) and is approximately 0.9% (right). 

is 0.1%. Both figures show that the modified NWA reduces the sizable deviations that 
occur in standard NWA as a daughter or parent mass is approached to the conventional 
uncertainty estimate. A multiple, overlapping application of Eq. (2.9) that would also 
include gluino production and chargino decay effects could be envisioned, but is beyond 
the scope of this work. 

4. Conclusions 

For configurations with kinematical bounds in the vicinity of resonances phase space sup- 
pression via /3-f actors can significantly distort the resonance Breit-Wigner, thus effecting 
an unexpectedly large NWA error. For affected configurations we proposed a modification 

2 The chargino is treated as stable and the gluino in standard NWA with spin correlations. 
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of the standard NWA that allows to take this kinematical phase space suppression into 
account and thus to reduce the approximation uncertainty to the inverse of the generic res- 
onant enhancement M/Y. For supersymmetric extensions of the SM we have demonstrated 
this uncertainty reduction for similar masses in processes with single particle or cascade 
decay. If applied in phenomenological studies and data analysis with tools like Fittino [12], 
SFITTER [13] or MARMOSET [14], the method would contribute to an accurate deter- 
mination of BSM model parameters and thus to establishing supersymmetry or other key 
properties of the fundamental theory. 
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